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The effects of halo alignment and shape on the clustering 
of galaxies 
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ABSTRACT 

We investigate the effects of halo shape and its alignment with larger scale structure 
on the galaxy correlation function. We base our analysis on the galaxy formation mod- 
els of Guo et al., run on the Millennium Simulations. We quantify the importance of 
these effects by randomizing the angular positions of satellite galaxies within haloes, 
either coherently or individually, while keeping the distance to their respective cen- 
tral galaxies fixed. We find that the effect of disrupting the alignment with larger 
scale structure is a ~ 2 per cent decrease in the galaxy correlation function around 
r pa 1.8/i -1 Mpc. We find that sphericalizing the ellipsoidal distributions of galaxies 
within haloes decreases the correlation function by up to 20 per cent for r < 1 h^ 1 Mpc 
and increases it slightly at somewhat larger radii. Similar results apply to power spec- 
tra and redshift-space correlation functions. Models based on the Halo Occupation 
Distribution, which place galaxies spherically within haloes according to a mean ra- 
dial profile, will therefore significantly underestimate the clustering on sub-Mpc scales. 
In addition, we find that halo assembly bias, in particular the dependence of clustering 
on halo shape, propagates to the clustering of galaxies. We predict that this aspect of 
assembly bias should be observable through the use of extensive group catalogues. 

Key words: cosmology: theory - cosmology: large-scale structure of Universe 
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1 INTRODUCTION 

Investigating how matter is organized in our Universe is one 
of the key ways in which we can test the validity of cosmolog- 
ical models and constrain their parameters. By comparing 
theoretical predictions to observed measures of structure, 
such as the galaxy correlation function or the matter power 
spectrum, one can reject some models and fine-tune others. 
It is, however, important to keep in mind the limitations 
of theoretical models, both numerical and analytical, when 
making this comparison, as these may limit the applicability 
of the results. 

There are various ways in which one can predict the 
organization, or "clustering", of matter and galaxies the- 
oretically. One can use fully hydrodynamical simulations, 
in which dark matter, gas and stars are treated explicitly, 
to follow the formation and evolution both of dark matter 
haloes and of the galaxies within them. For a recent review of 
the nu merical methods behind such simulations, see lSpringell 
(J2010T) . Such models are computationally expensive, limited 
to small volumes in comparison to recent galaxy surveys, 
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and sensitive to the ad hoc subgrid recipes required to in- 
clude critical processes like star formation and feedback. 

An alternat ive, first imple mented bv jKauffmann et al.l 

(J1999T) (see also ISpringel et al . 2001, 2005), is to combine 
N-body simulations of the growth of dark matter struc- 
tures with sem i- analyt ic models of galaxy formation (e.g. 



White fc Frenkl 
19931 . ICole et all 



1991 . iKauffmann. White fc Guiderdonil 



1994 



Baughl l2006f ~for 



a review 



A 

great advantage of semi-analytic simulations is that they 
require comparatively little CPU time even for a large 
underlying N-body simulation. This allows them to be run 
many times and on many haloes, so that one can explore the 
physical processes and the associated parameters that are 
required to produce galaxy populations in agreement with 
selected observational data (such as the galaxy stellar mass, 
luminosity or correlation functions). Such semi-analytic 
simulations do not focus on the properties of individual 
objects, but rather on the underlying statistical properties 
of the entire population. In this way, the relative impor- 
tance of different physical processes can be examined as a 
function of the time and place where they are occurring. 
A disadvantage of such simulations is that they provide 
only very crude information on the structure of individual 
objects. 
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Yet another alternative is to take the statistical ap- 
proach one step further. If one is interested only in the 
present-day clustering of galaxies, the physical processes as- 
sociated with their formation and evolution may not be rel- 
evant. One can then populate the haloes in an N-body sim- 
ulation with galaxies using a purely statistical model that 
depends on current halo properties, for example halo mass. 
Galaxy clustering can then be described in terms of the 
clustering of the haloes. This approach is known as halo 
occupation distribution modellin g, or simply HOD mod- 
elling (see ICoorav fe Shethl 12002 for a review) . Typically, 
central and satellite galaxies are treated separately, as each 
halo will contain one and only one of the former but may 
contain none or many of the latter (|Kauffmann et al.lll993 . 
iKravtsov et alj2004IZheng et alj2005J ). The satellite galax- 
ies assigned to a halo are usually assumed to be spheri- 
ca lly distributed following a stand ard profile such as that 
of iNavarro. Frenk. fc White] l|l997f ). Attempts at including 
substructu re or an environmental dependence have also been 
made (e.g. iGiocoli et al.ll2010l . iGil-Marin, Jimenez fc Verde! 
I2011T ). Note that by assuming spherical symmetry some in - 
formation is lost. As the simulations of iDavis et al.l (|1985T ) 
first showed, cold dark matter haloes are typically strongly 
ellipsoidal. If the distribution of galaxies follows the mass 
distribution, this would leave an imprint on the galaxy cor- 
relation function on small scales. Furthermore, halo elliptic- 
ity may also have an effect on larger scales. If neighbouring 
haloes are aligned, as expected from tidal-torque theory, this 
will boost the correlation on scales corresponding to the typ- 
ical separations between haloes. 

The ellipticity and intrinsic alignment of dark mat- 
ter haloes and their galaxy populatio ns have been the 



subjec t of many ear li er studies, e . g. 



ICarter fc Metcalfe 



J1980h.lBmggeiil dl982h . |Westl d 19891). ISplinter et all (|l997t ). 
IJing fc Sutoi J2002h an d bailin fc Steinmet? d2005l) . and re- 
centlv lPaz et al.l (|201ll ') and lSmargon et alJJ2012h. Most rel- 
evant to the current work ar e the studies bv lSmith fc Wattsl 
d2005l) and IZu et al.l (|2008h . The former authors investi- 
gated the effects of halo triaxiality and alignment on the 
matter power spectrum in the halo model framework. In- 
spired by the results of simulations, they took a purely ana- 
lytic approach in which they re-develope d the hal o mode l 
to account for ellipsoidal halo shapes. IZu et al.l |2008), 
on the other hand , used the semi-analytical models of 
|Pe Lucia fc BlaizotJ 120071 ) to investigate environmental ef- 
fects, including that of halo ellipticity, on the galaxy corre- 
lation function in both real and redshift space. 

There exists a deeper connection between halo shape 
and cluste r ing t hat we also explore in thi s pape r. As 
iBett et all (|2007ft and iFaltenbacher fc White! (|2010h have 
previously shown, at fixed mass the clustering of haloes de- 
pends on their shape. This is probably a reflection of as- 
sembly bias (i.e. the dep endence of halo clus tering on prop- 
erties other than mass, iGao fc White! 120071 ). More specif- 
ically, the most spherical haloes in their samples cluster 
significantly more strongly than average, and the most 
aspherical more weakly. The known correlations between 
formation time and halo sh ape (e.g. lAllgood et al.l 120061 . 
iRagone-Figueroa et al.l 120101 ) , and between formation time 
and clustering strength (e.g. |Gao. Spr ingcl & White! 1200 
Wechsler etail |2006|" IWetzel et al 1 120071. IJing. Suto fc M 



be strong enough to explain the magnitude of the effect. 
Here, we investigate whether this shape-dependence is also 
seen in the clustering of the galaxies. If so, this would bring 
us one step closer to measuring assembly bias directly in 
observat ions. W e are h erein also motivated by the result s 
from IZhu et all (|2006h and ICroton. Gao fc White! (|2007h . 
who showed that assembly bias in general is indeed expected 
to propagate to galaxy clustering. 

In this paper, we expand upon previous work by in- 
vestigating the effects of alignment and ellipticity on the 
galax y correlation functio n using the Millennium Simula- 
tion iS pringcl ct al. 2005) and the semi-analytic models of 



IGuo et al 



(|2011T ) . In Section [2] we discuss these simulations 
and our methods for quantifying the effects of alignment and 
ellipticity. We also outline our procedure for determining the 
shape-dependence of galaxy bias. We show our results in 
Section |31 and present our conclusions in Section 3] 



2 METHODS 

2.1 Simulation and SAM 

We make u se of the galaxy catalogues generated 
by IGuo et al.l (|201lL hereafter Gil), who implemented 
galaxy formation models on the Millennium S imulations 
dSpringel et~aill2005l . iBovlan-Kolchin et al.ll2009l ). The Mil- 
lennium Simulation (MS) is a very large cosmological N- 
body simulation in which 2160 3 particles were traced from 
redshift 127 to the present day in a periodic box of side 
500 hT 1 Mpc, comoving. The Millennium-II Simulation (MS- 
II) follows the same number of particles in a box of side 
100 ft -1 Mpc and so has 125 times better mass resolu- 
tion. Both simulations assume a ACDM cosmology with 
paramete rs based on a combined analysis of the 2dF- 
GRS dColless et al.l l200ll ) and the first-year WMAP data 



|Spergel et al.ll2003 ). These cosmological parameters, given 
by {Q m , fib, n A ,as,n s ,h} = {0.25, 0.045, 0.75, 0.9, 1.0, 0.73}, 
are not consistent with the latest analyses of the CMB data, 
for example the seven-year WMAP results dKomatsu et al.l 
2011). In particular, the more recent data prefer lower as 
and higher f2 m valuesjj We will only make relative compar- 
isons between clustering statistics here, and do not expect 
our results to be signifi cantly influenced by these small pa- 
rameter differences (see lGuo et al.ll2012l ). 

The galaxy formation models of Gil allow galaxies to 
grow at the potential minima of the evolving population of 
haloes and subhaloes in the simulations. Each Friends-of- 
Friends (FoF) group contains a central galaxy at the poten- 
tial minimum of its main subhalo, and may contain many 
satellite galaxies at the centres of surrounding subhaloes. 
In some cases, due to tidal effects, a satellite galaxy may 
be stripped of its dark matter to the point where its sub- 
halo is no longer identified as a bound substructure, turn- 
ing the galaxy into an "orphan". Such galaxies follow the 
orbit of the dark matter particle that had the highest bind- 
ing energy immediately before subhalo disruption, except 
that their distance to the central galaxy is artificially de- 
creased until they merge with it in order to mimic the ef- 
fects of dynamical friction. We note that the treatment of 



2007J) point in this direction but may by themselves not 



1 See lAngulo fc White! ||2010| ) for a method to correct for this. 
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the orbits of orphans is approximate, and for example does 
not include the expected circularizati on of the orbits (e.g. 
iBovlan-Kolchin. Ma fc Quataertll200gf ). The models also in- 
clude treatments of star formation, gas cooling, gas strip- 
ping, metal enrichment, supernova and AGN feedback, and 
galaxy mergers. For more details about the SAM, as well as 
the treatment of different types of galaxies, we refer to Gil. 
For our purposes, it is enough to note that the predicted 
clustering of galaxies is quite a close match to that seen in 
the Sloan Digital Sky Survey iJGuo et al.1120111 ). 



2.2 Calculation of the galaxy correlation function 

The galaxy two-point correlation function, £(r), measures 
the clustering of galaxies as a function of scale. It effectively 
encodes the excess probability of finding a pair of galaxies at 
a given separation r, relative to the expectation for a uniform 
random distribution. In what follows, we will be interested 
in scales 30 h~ kpc < r < 50/i^ 1 Mpc, as these are both 
well-resolved and well-sampled by the simulation. In order 
to get accurate results over this full range, we calculate the 
correlation function by direct pair counts on small scales 
(i.e. r < 4/i~ 1 Mpc) and use an approximate but accurate 
method to calculate it on intermediate and large scales. 

A direct calculation of this function scales as the num- 
ber of galaxies squared and is thus unfeasible for the large 
sample analysed here. We therefore speed up the calcu- 
lation by mapping galaxies onto a grid, and we correlate 
the mean density contrast in each grid cell with that of 
every other (a method previously employed by, for ex- 
ample, Barriga fc Ga ztanaga 20021. lEriksen et al.l 12004 and 
ISanchez. Baugh fc Angulal2008l ). We improved the perfor- 
mance on intermediate scales by folding the density field 
onto itself before i ts autocorrelation is calculated (see e.g. 
I Jenkins et al.lll998T ). We do not go into these methods here, 
but note that tests against higher-accuracy calculations 
show that the error in the ratio of the correlation functions, 
which is the relevant quantity for our main results, is less 
than 1 per cent on all scales considered. For the correlation 
functions we calculate to determine the galaxy bias a direct 
pair count over the full range of scales is feasible, as there we 
only consider relatively small subsets of galaxies (see SfjO} . 



2.3 Testing the importance of alignment and 
ellipticity 

During their lifetime, haloes merge and may accrete 
more subhaloes. The accretion of mass is not isotropic 
sinc e matter flows in preferentially along filaments (s ee 
e.g. iTormen, Bouchet fcWliitell 19971 . IColberg et aill 19991 or 
more recently I Vera-Ciro et al.ll201 if ). As a result, the distri- 
bution of subhaloes and thus galaxies within a FoF group is 
generally not isotropic either, but is instead approximately 
ellipsoidal, following the mass and aligning with surr ound- 
ing large-scale structure (see e.g. lAngulo et al,ll2009h . To 
test whether alignment with neighbouring structure has an 
effect on clustering statistics, we randomly rotate the haloes 
around their centres and see if this systematically alters the 
galaxy correlation function. More precisely, we rotate the 
satellite population of each FoF group bodily around the 
central galaxy to a new randomly chosen orientation, and 



we repeat this process for every FoF group in the simulation. 
We stress that this transformation preserves the numbers, 
properties, and relative positions of the galaxies in every 
halo; only the orientations of the distributions change. We 
then calculate the galaxy correlation function for the new 
distribution, and compare it to the original. If alignment 
with large-scale structure is important, one would expect to 
see the correlation decrease systematically on scales slightly 
larger than individual haloes. To estimate the uncertainty 
in our results, we have repeated this process 25 times, each 
time with a different set of randomly chosen angles. 

The effect of halo ellipticity is tested in a similar way. 
Here, we randomly rotate the position of each individual 
satellite galaxy around its central, rather than rotating 
all satellites together. In this way, the galaxy distribution 
within each halo is sphericalized. Since the distribution of 
galaxies within haloes is typically ellipsoidal, this process 
should increase the average distance between galaxies, thus 
decreasing the co rrelations betwe en galaxies in the same 
halo. We note that lZu et al.l (|2008l ) investigated the effect of 
halo ellipticity in the same way. 

2.4 Testing the dependence of galaxy bias on halo 
shape 

iFaltenbacher fc White! l|2010h showed that the clustering 
of haloes depends on the shape of the halo, defined as 
s — c/a, where c and a are eigenvalues of the inertia tensor 
(a > b > c). Specifically, they showed that the large-scale 
bias of haloes with more spherical shapes is larger than aver- 
age, while the inverse is true for the most aspherical haloes. 
They also found that this difference decreases with equiv- 
alent peak height, i/(M,z) = S c (z)/a(M,z), where a(M,z) 
is the root-mean-square linear overdensity within a sphere 
which contains the mass M in the mean, and 5 c (z) is the 
linear overdensity threshold for collapse at redshift z. Here, 
we are interested in seeing if this shape-dependent cluster- 
ing, which might reflect the assembly bias of the haloes, is 

also recovered from the galaxy distri bution. 

W e use the halo shape data from IFaltenbacher fc White! 
(2010), who calculated the inertia tensor from the dark mat- 
ter particles belonging to the most-pronounced subhalqj 
of each FoF halo, which on average comprises ~ 80 per 
cent of its mass. To ensure that the shapes were accurately 
determined, only haloes with at least 700 particles were 
considered, corresponding to a minimum (sub)halo mass 
M = 6.02 x 10 11 h' 1 Mq. We compare this to the shape 
measured from the galaxy distribution in the same way, us- 
ing all galaxies with stellar masses M» > 10 9 /i -1 Mq within 
a sphere of radius R200 , defined as the radius enclosing 200 
times the mean density of the Universe, centred on the cen- 
tral galaxy. Using such a distance cut makes it easier to 
compare our results to observations - in fact, a similar pro- 
cedure is often followed when determining the richness of 

2 The most-pronounced subhalo may differ from the most mas- 
sive subhalo only if the FoF groups hosts two or more subhaloes 
with roughly the same mass. In this case one of these is arbitrarily 
assigned to be the most massive. The most-pronounced subhalo 
is more consistently defined by using the halo merger tree. It is 
also the subhalo hosting the most luminous galaxy in the FoF 
group. 
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Figure 1. The effect of halo alignment on the galaxy correla- 
tion function. The x-axis shows the real-space separation r, while 
the y-axis shows the fractional difference between the correla- 
tion function after random bodily rotations are applied, and that 
of the original, unrotated sample. All simulated galaxies with 
Af* > 10 9 h _1 Mq from the z = catalogue of Gil have been 
used here. The bin size is roughly 0.07 dex. Each of the 25 thin, 
coloured lines represents a different set of random rotations, and 
the thick, black line shows the average of these. There is a clear 
signal around r ~ 1.8 /i —1 Mpc, where the correlation function is 
lowered by roughly 2 per cent. 




I 
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r [Mpc/h] 

Figure 2. The effect of halo ellipticity on the galaxy correlation 
function. The bin size, axes and lines are as in Figure[JJ although 
the bodily rotations have been replaced by random independent 
rotations of satellites around their central galaxies. A peak of 1 — 2 
per cent can be seen around r fa 3.5 ft -1 Mpc, but the largest 
effect is seen on small scales, where the correlation function is 
systematically lowered by up to ~ 20 per cent. Note also that 
the scatter has been greatly reduced relative to Figure [l] This is 
mainly due to the larger number of random rotations used when 
rotating satellites separately. 



real groups and clusters. Rejecting galaxies outside the virial 
radius slightly biases us to measure more spherical shapes, 
but we have checked that this effect is small and does not 
significantly affect our results. 

After splitting the galaxies by the shape of their 
halo (measured either from the dark matter or from the 
galaxies themselves) we determine the large-scale galaxy 
bias factor 6 Ea i for each subsample. Here, too, we follow 

1 J 71 

Faltcnbacher & White (2010), who in turn followed the ap- 



proach of iGao fc White! (120070 . The bias is computed as 
the relative normalization factor that minimizes the mean 
square of the difference log(£ gm ) — log(6 ga i£ mm ) for four bins 
spaced equally in logr in the range 6 < r < 20 h" 1 Mpc. 
Here £ mm is the dark matter autocorrelation function and 
£ gm is the cross-correlati on function of galaxies and d ark 
matter. Note that unlike iFaltenbacher fc White! (|2010T ) we 
are only interested in the results at z — 0. 



3 RESULTS 

3.1 Alignment and ellipticity 

We will first discuss our results for "bodily" rotations, 
which test the effect of halo alignment. Figure [l] shows 
the fractional difference between the correlation functions 
of the "rotated" and original samples, plotted against the 
real-space separation r. We have only used those galaxies 
from the catalogue generated by Gil that have a stellar 
mass M„ > 10 9 /i _1 Mq, as the Millennium Simulation is 
not complete below this limit. This provides a sample of 
5 200 801 galaxies. We note that increasing this mass limit 
by a factor of ten does not influence our results, either qual- 
itatively or quantitatively. All random rotations are applied 
prior to the mass cut in order to avoid problems in cases 



where a central galaxy below the limiting mass has satellites 
above it. Coloured lines indicate different sets of rotations, 
while the thick, black line shows the average of these. There 
is clearly a significant dip around r fa 2/t -1 Mpc, with a 
depth of 2 per cent. This is due to the disruption of the 
alignment between haloes and surrounding structure. Note 
that the scatter is extremely low, due to the large number 
of objects (in fact, the uncertainty at large scales is domi- 
nated by the errors due to our approximate calculation of 
the correlation). Neglecting the orientation of haloes when 
populating them with galaxies will therefore have a modest, 
but significant, effect on the derived correlation function. 

The ellipsoidal shape of the haloes, and the fact that 
the galaxy distribution follows this shape, is a more signifi- 
cant factor when modelling the galaxy distribution. Figure[2] 
shows the result of applying independent rotations, which 
sphericalize the galaxy distributions within haloes. This sub- 
stantially suppresses correlations for r < 2/i~ 1 Mpc, with 
a ~ 20 per cent effect on the smallest scales probed here. 
The ellipsoidal shape of the galaxy distribution within haloes 
significantly reduces the typical separations of pairs within 
them. This is compensated by a 1 — 2 per cent stronger cor- 
relation around r ~ 3.5 /i _1 Mpc. We conclude that models 
that assume spherical profiles for the distribution of galax- 
ies within haloes will underestimate the galaxy correlation 
function by up to ~ 20 per cent, depending on the smallest 
scale considered. Note that the scatter is even smaller than 
before on all scales, which is due to the increased number of 
degrees of fr eedom her e . This result is in excellent agreement 
with that of IZu et al.l (I2008|). who appl i ed sim ilar methods 
to all galaxies from lDe Lucia fc Blaizotl (|2007l ) with r-band 
luminosities M r < — 19. 

One might worry that the necessarily artificial treat- 
ment of orphan, or "type 2" , galaxies in the galaxy formation 
models of Gil influences these results. Gil already showed 
that the inclusion of the orphans is critical if the radial dis- 
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Figure 3. Same as Figures [T] and [2] but now for the fractional differences in the galaxy power spectrum versus the wave number k. Left: 
Result when testing for alignment. Again a weak but systematic signal of a few per cent can be seen, now between k m 0.1 /iMpc and 
k « 2hMpc -1 . Right: Result when testing for ellipticity. A monotonic decline in power sets in at k ss 0.1 /i Mpc -1 , reaching roughly 
20 per cent at fc = 30/iMpc , matching the result found in Figure [2] This again demonstrates the importance of taking the ellipsoidal 
shape of the galaxy distribution within haloes into account. 



tribution of galaxies within rich clusters in the Millennium 
Simulation is to agree both with observations and with the 
much higher-resolution MS-II. We have investigated effects 
on our analysis by repeating it with these galaxies removed, 
reducing our sample size by ~ 24 per cent and spoiling the 
relatively good agreement of its small-scale correlation with 
observation. This removal significantly amplifies the signal 
found for the effects of alignment. This is because the orphan 
galaxies are primarily located near halo centres. Once they 
are removed, galaxies that do contribute to the alignment 
signal receive more weight. Orphan removal also changes 
the signal found for the effects of ellipticity, modestly boost- 
ing it down to r ~ 0.1 ft -1 Mpc. As a further check that 
the distribution of orphans in the simulation is realistic, 
we have examined the shapes of the galaxy distributions of 
massive haloes (specifically, 14 < log 10 (MF O F/[/i -1 Mq]) < 
14.5) in the MS and MS-II, again using Gil's galaxy cat- 
alogues and considering only galaxies with stellar masses 
M, > 10 9 h' 1 Mq. The better mass resolution of the MS-II 
results in far fewer orphans in this mass range, and conse- 
quently the positions of galaxies in MS-II are determined 
more accurately. Nevertheless, the shapes of the galaxy dis- 
tributions agree very well, thus implying that the distribu- 
tion of orphans in the MS is consistent with the distribu- 
tion of similar, but unstripped, galaxies in the MS-II. We 
also found that these shapes agree ver y well with thos e of 
the dark matter haloes themselves fsee lBett et al.ll2007l ). A 
more detailed discussion of the shapes is beyond the scope 
of this paper. 

For completeness, we have also compared the galaxy- 
galaxy power spectra of the rotated and unrotated samples. 
The results are shown in Figure [3] Three foldings in to- 
tal were used to calculate the power spectra over the full 
range shown, each with a fold factor of six (i.e. each folding 
maps the particle distribution to l/216th of the volume). 
The power spectra were re-binned logarithmically to resem- 
ble the bins used for the galaxy correlation functions, and 
to reduce noise. Shot noise, which dominates the power for 
A; > 10 /i Mpc -1 , was subtracted. The left-hand figure shows 
the fractional differences that result from applying bodily ro- 



tations. Just as for the correlation function, there is a weak 
but clear dip of 1 — 2 per cent present which reflects the align- 
ment of haloes with surrounding structure. As expected, the 
right-hand figure shows a much stronger decrease in power, 
up to 20 per cent on the smallest scales considered. This 
again demonstrates the importance of taking the ellipsoidal 
distribution of galaxies within FoF groups into account in, 
for example, models that use that use the full shape of the 

power spectrum to extract cosmological pa rameters. 

Our results differ from those found bv lSmith fc Watts! 

(2005). Like ourselves, they find that the scale at which the 
contribution from alignments to the power spectrum is max- 
imal is ~ 0.5 h Mpc -1 , but they show that the relative con- 
tribution of alignments is strongly model-dependent, varying 
from 10 -12 to 10 per cent. Furthermore, they find that when 
haloes are assumed to be spherical, the power is higher than 
when they are ellipsoidal by up to 5 per cent. Not only is 
the effect we find significantly stronger and increasing to- 
wards smaller scales up to at least k = 30 ft Mpc -1 , but its 
sign is opposite. We attribute t hese differences to the fact 
that in both models explored bv lSmith fc Watts! i|2005r ), the 
radially averaged density profiles are not conserved when 
transforming the haloes from spherical to triaxial, making a 
comparison with our own results difficult. 



3.2 Shape-dependent galaxy bias 

Having established that the shape of the galaxy distribution 
significantly affects the small-scale clustering, we now inves- 
tigate how the shape-dependent assembly bias affects clus- 
tering. We first examine how well the shape measured from 
the galaxy distribution corresponds to that measured from 
the dark matter particles belonging to the most-pronounced 
substructure. Note that while the galaxies trace the shape 
of the FoF halo very well, this is not necessarily the case 
for its most-pronounced substructure. Additionally, even if 
the galaxy distribution traces the dark matter perfectly, the 
right shape may still not be recovered if only a small number 
of galaxies is available to sample the halo. 

The results of this shape comparison are presented in 
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Figure 4. A comparison of the halo shape measured from the 
dark matter and that measured from the galaxy distribution. Here 
shape is defined as the ratio of the smallest to the largest eigen- 
value of the inertia tensor. The shape on the horizontal axis is 
computed using all galaxies satisfying M* > 10 9 h _1 Mq within 
R200- For the shape on the vertical axis, only CDM particles be- 
longing to the most- pronounced substructure are included, and 
only if the number of particles is at least 700. Pixels are colour- 
coded by mean number of galaxies per halo. The "true" shape is 
recovered more accurately when the number of galaxies increases; 
some scatter remains however, as the galaxy distribution does not 
perfectly trace the main subhalo. 




i/(M) 
Figure 5. The dependence of galaxy bias, b ga i, on shape, as a 
function of peak height at z = 0. The black line in the top panel 
shows the galaxy bias of the full galaxy sample, while coloured 
lines show the bias of subsamples split by halo shape measured 
from the dark matter distribution. In the bottom panel the frac- 
tional differences of the bias of these subsamples relative to the 
full sample are shown. The vertical error bars show ltr devia- 
tions calculated from 50 bootstrap resamplings of the galaxy cat- 
alogues. We find that galaxies in the most spherical (aspherical) 
haloes are strongly biased (antibiased) relative to the full sam- 
ple. The difference can be as much as 40 per cent for v K, 0.7. 
For v > 2.5 statistical uncertainties, due to the low number of 
high-mass haloes, begin to play an important role. 



Figure [4] Here the shape measured from the galaxy distri- 
bution is shown on the horizontal axis, while the shape mea- 
sured from the dark matter is on the vertical axis. Each pixel 
is colour-coded by the mean number of satellites, iV sa t, sat- 
isfying M» < 10 9 /i _1 Mq and |r sat — r ccn | < i?200- A low 
value of s = c/a indicates that the halo is (measured to 
be) very aspherical, while a perfectly spherical halo would 
have s — 1. It is immediately clear that a low number of 
satellites leads to a severe underestimate of s. This is ex- 
pected: together with the central galaxy, any two satellites 
will define a plane, ensuring that c = 0. It is only when the 
halo is sampled by a large enough ensemble of points that 
the inertia tensor can be determined accurately. Figure [4] 
illustrates that one needs 7V sat > 30 to get an unbiased and 
accurate shape estimate. However, there is always a signifi- 
cant amount of scatter around the diagonal. This is due to 
the galaxies tracing the mass, i.e. the shape of the whole 
FoF group, and not just that of the most-pronounced sub- 
structure. 

Next, we split our galaxy sample by the shape mea- 
sured from the halo dark matter distribution and calculate 
the galaxy bias factor as function of equivalent peak height, 
fegai(^), following the method described in £|2.4I As we can 
only use galaxies for which a dark matter shape has been 
determined - i.e. those in FoF haloes of which the most 
massive substructure is comprised of at least 700 particles - 
our galaxy sample is reduced to 2 953 050 galaxies. The re- 
sults are shown in Figure [5] Here the upper panel shows the 
galaxy bias determined for each sample at a given equivalent 
peak height. The points show the median value, horizontal 



error bars indicate the width of the bin, and vertical error 
bars show la deviations calculated from 50 bootstrap re- 
samplings of the galaxy catalogues. The black line shows 
b S ai(v) for the full sample of galaxies, while the coloured 
lines show the bias for the different subsamples. In the bot- 
tom panel, the fractional difference between these subsam- 
ples and the full sample is shown. It is immediately clear that 
there is a strong dependence on shape: at the lowest equiv- 
alent peak heights probed here, the bias of the galaxies in 
the most spherical (aspherical) 20 per cent of haloes is up to 
40 per cent higher (lower) than that of the full sample. This 
shows that the sh ape- dependence of the halo bias fo und by 
iBett et all (|2007h and IFaltenbacher fc White! IJ2010l l is also 
strongly present in the clustering of the galaxies. The effect 
grows weaker with increasing v. However, for v > 2.5 statis- 
tical uncertainties begin to dominate, due to the low number 
of haloes available at high equ ivalent peak heights. No t e that 
this is not seen in the data of IFaltenbacher fc White! (|2010f ) 
as they combine the data from different redshifts, while we 
only consider z = 0. 

We then repeat this exercise, but this time we split our 
galaxy sample by the shape measured from the distribu- 
tion of the galaxies themselves. However, we expect low- 
mass haloes to host only a few satellite galaxies, leading to 
very unreliable estimates of the halo shape (see Figure [4}. 
In order to separate the signal we are looking for - i.e. the 
shape-dependence of the galaxy bias - from the unwanted 
bias introduced by using too few galaxies in the shape mea- 
surement, we now consider fe ga i as a function of the number 
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Figure 6. As Figure [5] but now showing the bias as a function 
of the number of satellite galaxies and split by the halo shape 
measured from galaxies within ifeoo- At low iV sat no significant 
shape-dependence of 6 ga i is recovered, due to the extremely unre- 
liable shape determinations that follow from using only a handful 
of galaxies to sample the halo. For 10 < A^at < 400, however, 
we find again that galaxies in more spherical (aspherical) haloes 
have a significantly higher (lower) bias than average. At higher 
Ns&t our results are again dominated by poor statistics. 



of satellites per halo, iV sa t, instead of the equivalent peak 
height v. An additional advantage of this approach is that 
Af Ba t is directly observable. We note, however, that we ob- 
tain almost identical results when considering galaxy bias as 
a function of peak height instead of iV S at • Since we do not use 
the dark matter particle data in this case, we are no longer 
constrained by needing haloes with at least 700 particles. 
However, as we can now only consider haloes with at least 
one satellite galaxy that satisfies both M t > 10 /i _1 Mq 
and jr sa t — r ccn | < -R200, we are left with a sample of 2 566 441 
galaxies. The results are shown in Figure [6] At the lowest 
value of A^sat, no significant effect can be seen. But as A^at 
grows, increasing the accuracy of the shape determinations, 
we again see a clear dependence of 6 ga i on the shape s: galax- 
ies in more spherical haloes are ~ 20 per cent more strongly 
clustered than average. The inverse is true for galaxies in 
the most aspherical haloes. When A^at grows too high our 
results are once more dominated by statistical errors, due to 
the low number of high-mass haloes (hosting at least several 
hundreds of satellites) available. 

These results show that halo assembly bias in the form 
of a shape-dependent clustering strength propagates to the 
clustering of galaxies, and can therefore in principle be mea- 
sured in sufficiently large surveys. In order to carry out such 
a task, a large galaxy survey with appropriately defined 
group catalogues is needed. 



jSpringel et al.l 120051 , iBovlan-Kolchin et all 120091 ) and the 
galaxy formation models of iGuo et al.l 1 20111) . By rotating 
satellite galaxies in FoF groups around their central galax- 
ies, either coherently for each halo or independently for each 
satellite, and then comparing the correlation function of the 
resulting galaxy distribution to the original one, we were 
able to quantify the importance of taking halo alignment and 
non-sphericity into account. Furthermore, by measuring the 
shape of the haloes as traced by the galaxies we were able 
to investigate the propagation of shape-dependent assembly 
bias to the clustering of galaxies. Only galaxies with stellar 
masses M, > 10 9 h~ x Mq were considered in our analysis, 
though we note that increasing this mass limit by a factor 
of ten does not influence our results. Our findings can be 
summarized as follows: 

• The effects on the galaxy correlation function of the 
alignment of haloes with larger-scale structure are small. 
The main effect of disrupting this alignment is a 2 per cent 
reduction in correlation amplitude around r ~ 1.8 h~ x Mpc, 
with minor effects of at most 1 per cent at smaller scales. 

• The ellipsoidal shapes of the galaxy distributions within 
individual haloes have a much stronger influence on galaxy 
correlations. By sphericalizing these galaxy distributions 
(i.e. randomizing the angular positions of satellites while 
keeping the distance from the central galaxy fixed), the 
correlation function is raised by up to 2 per cent around 
r ~ 3.5 h" 1 Mpc, but greatly reduced for r < 1.5 h^ 1 Mpc, 
by up to ~ 20 per cent on the smalle st scale probed , 
r = 30 h' 1 kpc. This confirms the results of lZu et al.l l|2008h . 
The effect on the galaxy power spectrum extends to scales 
as large as k — 0.1/tMpc -1 . 

• The assembly bias of haloes, as characterized by the 
dependence of clustering on halo shape, is reflected in the 
clustering of galaxies. The effect is strongest at low equiva- 
lent peak heights: at v ~ 0.7, the galaxy bias of galaxies in 
the 20 per cent most spherical and most aspherical haloes 
deviate from the average by 40 per cent. 

• Even if the shape of the halo cannot be measured di- 
rectly, but is instead estimated from galaxies within one 
virial radius of the central galaxy, the effect of assembly 
bias is clearly visible. 

By using the plane-parallel approximation and ignoring evo- 
lution, we have checked that comparable results are obtained 
for the effects of alignment and ellipticity on the redshift- 
space correlation functions. Models that assume a spher- 
ically symmetric profile for the galaxy distribution, such 
as HOD models, will therefore significantly underestimate 
galaxy correlations and power spectra on sub-Mpc scales. 

Furthermore, we have demonstrated that the shapes of 
haloes and of the galaxy distributions within them can be 
strongly correlated with their clustering. This effect should 
be measurable in galaxy redshift surveys. With the help of 
extensive group catalogues it should therefore be possible to 
measure assembly bias directly. 



4 SUMMARY 

We have investigated the effects of halo alignment with 
larger scale structure and of halo ellipticity on galaxy 
correlation functions, using the Millennium Simulations 
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